Effective quantum dimer model for trimerized kagome antiferromagnet 
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An effective spin-orbit Hamiltonian is derived for a spin- 1/2 trimerized kagome antiferromagnet 
in the second-order of perturbation theory in the ratio of two coupling constants. Low-energy singlet 
states of the obtained model are mapped to a quantum dimer model on a triangular lattice. The 
quantum dimer model is dominated by dimer resonances on a few shortest loops of the triangular 
lattice. Characteristic energy scale for the dimer model constitutes only a small fraction of the 
weaker exchange coupling constant. 
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I. INTRODUCTION 

Resonating valence bond (RVB) statei is nowadays a 
popular paradigm in condensed matter physics. Short- 
range RVB states are considered to be probable candi- 
dates for an elusive spin-liquid phase of magnetic insu- 
lators. The idea of short-range RVB states is quanti- 
tatively formulated by so called quantum dimer mod- 
els (QDM)i2i&4 In a QDM each dimer represents a sin- 
glet state (valence bond) between a pair of neighboring 
spins. The QD Hamiltonian is defined in the Hilbert 
space of close-packed dimer coverings of a lattice. The 
dimer states are assumed to be properly orthogonalized. 
Local dynamics of an RVB state is typically described 
on the smallest plaquettes (□), which are squares for a 
square lattice or rhombi for a triangular lattice: 

^qd = £[-*(t>o H+H i)(=l)+^(l=> (=I+H i> (i il) 

□ 

... w 

The first term is a dimer kinetic energy, which flips a 
pair of parallel dimers around an arbitrary plaquette; the 
second term is a potential energy between such pairs. 

Rokhsar and Kivelson 3 - have shown that a short-range 
RVB state given by a superposition of all dimer cov- 
erings of a square lattice is an exact eigenstate of the 
QD Hamiltonian for a special choice of the parameters 
t = V. On a bipartite square lattice, the RVB state at 
the Rokhsar-Kivelson (RK) point has long-range power 
low correlations and describes, consequently, a gapless 
spin-liquid stateA^ Small perturbations away from the 
RK point drive the system into one of the ordered crys- 
talline dimer states. The QDM on a triangular lattice ex- 
hibits quite a different behavior at the RK pointi£*£ The 
short-range RVB state has exponentially decaying corre- 
lators and is fully gapped. It exists, therefore, in a finite 
range of parameters around the RK point and is stable 
with respect to weak perturbations to the QD Hamilto- 
nian. Still, question whether such states or Hamiltonians 
can describe realistic quantum spin systems remains un- 
settled. In the present work we propose a realization 
of QDM on a triangular lattice for a nearest-neighbor 
Heisenberg spin model. 

The most probable candidates for a singlet spin-liquid 




FIG. 1: Trimerized kagome lattice with two exchange con- 
stants: Ji in A-triangles and ,li in v _ triangles. 



ground state are frustrated quantum antiferromagnets. 8 
Numerical exact diagonalization studies of a spin-1/2 
Heisenberg kagome antiferromagnet have shown that this 
spin model has a nonmagnetic ground state with a large 
number of low-lying singlet excitationsAiS Accessible 
cluster sizes do not allow to draw a definite conclusion on 
the possible nature of the magnetically disordered (sin- 
glet) ground state. 

One of a very few analytic approaches to such strongly- 
correlated spin systems is an expansion from small clus- 
ters. The main motif of a kagome lattice is triangle. It 
is, therefore, natural to start from a trimerized kagome 
lattice shown in Fig.^ Such a strong-coupling approach 
has been pursued in relation to kagome antiferromag- 
net in several theoretical works Recently, an ex- 
perimental scheme to create a trimerized kagome lattice 
has been suggested for ultracold atomic gases in optical 
traps^ This opens a way for an experimental probe of 
RVB physics in the corresponding spin model. 

The Heisenberg model on a trimerized lattice 
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has two coupling constants: J\ for a stronger interaction 
between spins in A-triangles and Ji for a weaker interac- 
tion inside y-triangles. An array of isolated A-blocks is a 
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zeroth order Hamiltonian, which has a highly degenerate 
ground state. Interblock interaction lifts such a degener- 
acy. In section II we derive an effective Hamiltonian up to 
the second-order in a small parameter J = J2/J1 <C 1. 
This Hamiltonian is mapped to a QDM in section III. 
The obtained results and their implication for the ground 
state properties of the trimerized kagome model are dis- 
cussed in section IV. 



II. STRONG-COUPLING EXPANSION 

A. First-order Hamiltonian 

Let us in the beginning rederive the previous results on 
the effective first-order Hamiltonianii*i2ii^ using some- 
what different notations. Below we normalize all energies 
to J\ such that J 2 —* J~ In a strong-coupling expansion 
one starts with an isolated triangle described by 



Tl/\ — Si • S2 + S2 • S3 + S3 • Si 
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The energy levels of Ha are determined by the total 
spin 5 to t- For on-site S = 1/2, which is always as- 
sumed below, the levels are two degenerate doublets with 
5tot = 1/2 and E = — § and one quartet with Stat = 3/2 



and E 



The doublet states with 5f ot = 1/2 are 



dr) = -^(tU-UT), Ipt) = 77=(2m-TU-UT) 
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where spin numbering in an individual triangle follows 
Fig. [2k. The former state is a combination of the 
spin-up apex spin and a singlet bond between the two 
base spins and has odd parity under the permutation 
p23- The other state \p-[) is even under P23. The two 
members of a quartet with 5* ot = +3 and +1 are 

l?+s) = I TTT> , k+i) = -^(ITT + TIT + TTI) • (5) 

All other states are obtained from (J3J and (JSJ by applying 
5 t ~ t operator. 

The choice of the basis in the doublet subspace is not 
unique. The apex spin can be put into a singlet state 
either with its right or left neighbor. The two alternative 
basis states obtained by rotating |dj) about a center of 
triangle counterclockwise on 2ir/3 and 4ir/3 are 

14 



-^(ITT-TTI) , I4') = ^|(T1T- ITT) (6) 



with their orthogonal partners \pL) and \p"). Transfor- 
mation from the old basis J3J to the new states is 
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FIG. 2: Three different geometries of A-blocks contributing to 
the second-order energy correction in the interblock coupling. 
The labeling of axes and sites inside triangles is shown in the 
upper panel (a). 



where a =T, | or 1,2 is a spinor index. The main differ- 
ence with the previous worksii*i£ is that real basis states 
Q or J7J) are used instead of complex chiral states. This 
yields a more transparent form of the effective Hamilto- 
nian and simplifies subsequent derivation of a QDM. 

At this point we introduce two sets of the Pauli matri- 
ces: a 1 , which act between spin-up and spin-down states, 
and t 1 , which act in the orbital subspace (d,p) preserving 
the total spin. A convenient choice of orbital axes shown 
in Fig. corresponds to 



l K) = \d a ), r Zl \p a ) = -\ Pa ) 



(8) 



Then, the orbital operators projected onto the rotated 
axes 
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have the following eigenstates: 
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In order to find the effect of interblock coupling in 
the first order of perturbation theory in J = J2/J1 one 
should neglect S to t = 3/2 states separated by a finite gap 
AE = I and calculate matrix elements of the on-site spin 
operators between the low-energy doublet states \d a ) and 
\p a ). This problem is greatly simplified once all sym- 
metries are taken into account. Introducing operators 
dJjO) = \d a ) and pJjO) = \p a ), where |0) is a fictitious 
vacuum, the Hubbard-type representation of on-site spins 
is written as 

s i = 5 ^"^d - g PlvaPPp , (11) 
S 2: 3 = ^ pl<r a ppp ± —/= (pto-apdp + h. c.) . 

Spinor structure is a consequence of the spin-rotation 
symmetry, while permutation of the base spins P23 fixes 
the orbital part in (|1I|I . The above representation is fur- 
ther simplified once the total spin of a triangle S = \<j 
is defined and the orbital operators r Zk are used. Then, 



the 



,th 



spin (n = 1-3) of the i triangular block is rep- 



resented by 



^ni — g Sj (1 



27f) 



(12) 



where z n goes from the center of a triangle in the direc- 
tion of the corresponding spin, see Fig. 

The effective first-order Hamiltonian in the interblock 
coupling is found by substituting Eq. l(T2"|) into the Hamil- 
tonian (0): 



^ = : crE s ^ s ^ 1 + 2r ^)( 1 + 2r / m )' 



(13) 



where a trivial constant term — 4-/Va is omitted for con- 
venience. The derived spin-orbital Hamiltonian TL\ is 
defined on a triangular lattice, such that every site cor- 
responds to one A-block of the trimerized kagome model 
and is attributed with spin and orbital operators. The 
bond orientation uniquely determines the orbital axes for 
two participating sites. 



B. Second-order Hamiltonian 

The second-order energy corrections for weakly- 
coupled spin triangles have been obtained by Raghu 
and co-workers^ These authors have mostly been inter- 
ested in a one-dimensional model, therefore, their anal- 
ysis misses several terms relevant for a two-dimensional 
array of triangles in the trimerized kagome model. In or- 
der to calculate the second-order result in J 2 one needs 
to determine matrix elements of on-site spins between 
doublet (@J and quartet (J5J states. Introducing symmet- 
ric third-rank spinor tensor g tt| a 7 such that qm = 1 5+3) , 
9ii2 = 9i22 = and 9222 = \q-s) : and 



utilizing spin-rotation symmetry we find by analogy with 
Eq. GB 



Si =-^vUr{™ v ) fh P« + - ( 14 ) 
S2.3 = -^=ql Pl ((Ta y ) pry (±^-d a + ip a ) +h.c. 

The second-order energy correction in the interblock 
coupling is, generally, given by 



H 2 (G,G') = J2 



(G\H\X)(X\H\G') 



x 



En — E 



X 



(15) 



where G and G' denote combinations of lowest doublet 
states on A-blocks and X are excited states. The nonzero 
second-order terms appear if either (i) one J2-bond acts 
twice in the numerator of Eq. I|15|) or (ii) two adjacent 
J2-bond emerging from the same A-block are used subse- 
quently in the matrix elements (G\H\X) and (X\TL\G'). 
This determines three different geometries for two- and 
three-block interaction terms shown in Fig. [21 In the 
first case of two-block interaction, Fig. [2^, either one or 
both triangular blocks have quartets in the intermediate 
states X. For the three block interactions (Fig. Hb,c), 
only a middle block has excited quartets in the inter- 
mediate states. Every pair of free (uncoupled) spins in 
a A-block imposes an extra permutation symmetry on 
the second-order Hamiltonian H2 (G,G'). For example, 
the two-block cluster in Fig. [2K has extra Pi2 symme- 
try for the left z-block and P13 symmetry for the right 
/c-block. Therefore, the orbital states, \d'^) or \p'^) for 
the left block and \d' a ) or \p' a ) for the right block, remain 
unchanged during the second-order perturbation process 
l(ISfr. In other words 7i 2a commutes with r/ 3 and r^ 2 . 
The conservation of orbital state is also fulfilled for all 
triangles in the three-block term in Fig. |2t> and for the 
side triangles in Fig. [5J;, whereas the middle block does 
change its orbital state during the second-order process. 
The above conservation laws significantly simplify sum- 
mation over intermediate states in Eq. (|15|) . The final 
results are 
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(16) 



for the two-block interaction; 
4 7 2 

H *> = -243 E Si-S*(l+27? )(l+2Tf")(l-7f ) (17) 

(ijk) 

for the three-block interaction shown in Fig. [2t> with cor- 
responding labeling of blocks; and 

9 T 2 r 
Vic = — X)(l + 2r/-)(l + 2^)[s < .S k (l + 27f) 

(ijk) 



+ V3tV S 4 • (Sj x S fc ) 



(18) 
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for the three-block interaction in geometry of Fig. [2k- 
Polarization of orbital operators t?" in the above equa- 
tions is again found by simple inspection of the arrange- 
ment of corresponding blocks on a trimerized kagome lat- 
tice. Three-block terms in H.2b exist on y-plaquettes 
of an effective triangular lattice with three different z- 
axes in Eq. 1)17(1. Three-body terms in Ti.2 C appear on 
A-plaquettes of the triangular lattice with one polariza- 
tion of t z operators in Eq. I|18(l . which changes under 
permutation of (ijk). 

The interactions Jiia and H.2c coincide up to a triv- 
ial change of notations with the previously derived terms 
for a one-dimensional array of triangles^ while the term 
li.2b is a novel one. Note, that 7i.2 C contains a remarkable 
three-body spin-chiral interaction term. By deriving the 
effective spin-orbital Hamiltonians TL\ and TL2, we have 
substantially restricted the Hilbert space and simplified 
the original problem of finding the ground and the low- 
est energy states of the spin model Q. The remain- 
ing problem of solving Ti\ + Ti.2 is still highly nontrivial. 
The spin-orbital Hamiltonian Ti\ has been studied so far 
in the mean-field approach^ An effective Hamiltonian 
derived by the contractor renormalization group, which 
partially resembles Ti\ + Ti.2, has also been analyzed in 
the mean- field approximation. 14 Below we discuss a map- 
ping of the derived Hamiltonians lfH5|) and l(TH|l - (fT%|l to an 
effective QDM. The obtained QD Hamiltonian is domi- 
nated by the kinetic energy for dimer tunneling. The 
mean-field approximation, which assumes a frozen pat- 
tern of dimers, is, therefore, a poor approximation in the 
present problem. 



III. QUANTUM DIMER MODEL 

A. General remarks 

Search for the low-energy states of the first-order 
Hamiltonian i|13|l can be started by considering first a 
two-site problem (two adjacent A-blocks of the original 
kagome lattice)^ This problem is solved exactly and its 
ground state corresponds to a spin singlet with the orbital 
degrees fully polarized along the bond: (rf n ) = (tJ" 1 ) = 
1. The ground-state energy is — | J. A variational singlet 
ground state for the lattice problem H13[) is constructed 
by splitting the whole lattice into a close-packed struc- 
ture of dimers between nearest-neighbor sites, such that 
the dimer wave-function is given by the ground-state so- 
lution of the two-site problem. A remarkable feature of 
these variational states is that at the mean-field level with 
respect to orbital degrees of freedom the total energy is 
just a sum of energies of individual dimers and does not 
depend on a particular dimer covering of a triangular 
latticed Indeed, once (r t ? n ) = 1, then for the two other 
axes (r* m ) = (rf*) = —k. Therefore, the expectation 
value of any empty bond, i.e., a bond without a dimer, 
identically vanishes over the variational wave-function: 




FIG. 3: Effective triangular lattice with five shortest loops. 
The arrow directions indicate the sign convention for singlet 
wave-functions on each bond 

either one or both sites of the bond have (1 + 2r? m ) = 0. 

The degenerate set of variational mean-field states 
has been identified with low energy states of spin-1/2 
antiferromagnets on trimerized and isotropic kagome 
latticesii2*i& The number of low-lying singlets of the 
kagome model scales, then, as 1.15^ in good agreement 
with the full exact diagonalization study^fi The previous 
works leave, however, without answer question about va- 
lidity of the mean-field approximation and further lifting 
of degeneracy by quantum fluctuations. 

In order to beyond the mean-field approximation, 
one has to consider off-diagonal matrix elements of the 
Hamiltonian l|13|l between various dimer configurations 
as well as the corresponding overlap matrix. The general 
rule to compute the overlap matrix for models, where 
every dimer represents a singlet pair, is to construct 
transition or overlap graph by drawing two dimer con- 
figurations on the same latticed Every closed noninter- 
secting loop of dimers contributes 2/2'/ 2 to the over- 
lap matrix, I being the length of the loop. The sign of 
the overlap matrix element depends on a sign conven- 
tion for singlet wave- functions. We adopt the standard 
convention 2 *^ such that the singlet bond wave-function 
is [ab] — -^(lali, — latb), where b is an upper site in the 
pair or is directly to the right from a, see Fig. |3J 

Local dynamics of singlet bonds in trimerized kagome 
model is determined by a few shortest loops on an ef- 
fective triangular lattice, which include two- and three- 
dimer moves, see Fig. [3] Taking into account the orbital 
part of the wave- functions the overlaps of two dimer con- 
figurations on each loop are calculated as c\ = — 1/2 4 , 
c 2 = c 3 = -1/2 7 , c 4 = -1/2 8 , and c 5 = 1/2 5 . These 
overlap matrix elements are significantly smaller than 
for singlet bond configurations on the original triangu- 
lar lattice. In the latter case the corresponding loops 
have ci = 1/2, c 2 = c 3 = 1/2 2 , c 4 = c 5 = -1/2 2 . The 
difference reflects the fact that loops on an effective tri- 
angular lattice correspond to significantly longer loops on 
the original trimerized kagome lattice. For example, the 
shortest C\ loop corresponds to a loop of length I = 10 
on a kagome lattice. Loops C4 and C5 are different for 
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the considered model because kagome lattice has only a 
three-fold rotation axis in the center of every triangle. 
Significant difference of the overlap matrix elements ex- 
plains why a QDM description is a poor approximation 
for a spin-1/2 Heisenberg antiferromagnet on a triangu- 
lar lattice^ but may be a good one for the trimerized 
kagome antiferromagnet. Below in this section we assume 
that the ground states of the first and second-order effec- 
tive Hamiltonians are given variationally by close-packed 
dimer states and compute a new effective QD Hamilto- 
nian. The above assumption is supported by numerical 
treatment of the trimerized kagome antiferromagnet^ 

Derivation of a QDM from a particular spin Hamilto- 
nian has been formulated via calculation of the inverse 
square root of the overlap matrix^ We find that actual 
calculations become more transparent by operating with 
the wave-functions. The final result are, of course, equiv- 
alent in both approaches. Specifically, let us consider two 
linearly independent, normalized states \ipi} and \if)2), 
which have a small overlap (tpi\ip2) — (^IV'i) — c - Ma- 
trix elements of the Hamiltonian between the two states 
are assumed to be known En = (ipi\H |t/>i), E 2 i = E i2 = 
(-01 1 H\tp2), and E22 = (ip2\H\ip2)- The aim is to compute 
matrix elements in a new properly orthogonalized basis 
\ip n ). Transformation to the new basis is given by a sym- 
metric matrix: 

|y>i) =A(|Vi)-/#2», \V2) = A(IV*> - ■ (19) 

Conditions (ipi\f2) = 0, (<fi\tpi) = {^2^2) = 1 deter- 
mine /i and A. Assuming En = E22, for simplicity, and 
calculating matrix elements of 7i between the new states 
one finds 

En = En + - ^(cEix — E12) , 

1 — c/ 

E12 = E V2 + —^(cEi2 - En) • (20) 

In the following we shall subtract from the effective 
Hamiltonians Tli and Tii + TL2 the corresponding mean- 
field energies £mf, which are the same for all dimer cov- 
erings. Then, En = E22 = and the matrix elements 
(120(1 are directly related to the parameters of a QDM: 

t = -E12 « -E12 , V = En « -cE u , (21) 

where |c| <C 1 is used. 

B. First-order mapping 

Let us apply the outlined procedure to the first-order 
Hamiltonian l(13|) . The mean-field energy Emf = —jJNd 
of an arbitrary configuration of N e i dimers is always sub- 
tracted from Eq. 1)13(1 . The wave-functions for two dimer 
states on the shortest loop Ci shown in Fig. El are written 
explicitly as 

= [12] [43] |4444>, l^2> = [32] [41] 14444), 

(22) 



where sites are numbered counter-clockwise beginning 
with the lower right corner of the rhombus. The first part 
of 1^1,2) is given by a product of two spin singlet states, 
whereas the second part is an orbital wave- function repre- 
sented as a product of states (QJ or J7J. Explicit calcula- 
tion of the off-diagonal matrix element E12 — (V>2|7^i|0i) 
yields 

h = -^J , Vi = ^ J . (23) 

The largest kinetic energy term in the QD Hamiltonian 
amounts to less than 5% of the weaker coupling con- 
stant. The ratio of the potential energy to the kinetic 
term constant is also very small V/\t\ = 1/16. Note, 
that t < from the above calculation. Remaining free- 
dom in the choice of sign of tunneling matrix elements is 
discussed in the next section. 

Since the potential energy Vi is an order of magnitude 
smaller than the dimer hopping t±, the next relevant in- 
teractions besides the kinetic energy of two-dimer moves 
around Ci may be three-dimer resonances along longer 
loops C2, C3, C4, and C5. We find no tunneling for loops 
C4 and C5, i.e., V^s = £4,5 = 0. For loop C5, vanish- 
ing of the off-diagonal matrix element can be understood 
by drawing two dimer configurations on a corresponding 
cluster of a kagome lattice, which is a six-point star. Two 
dimer states around perimeter of such a star are exact de- 
generate eigenstatesii and, hence, t$ = 0. In the former 
case, loop C4, the E12 — result is valid only in the first 
order in J, see the next subsection. 

Coherent motion of three dimers along composite loops 
C 2 and C 3 is not described by Eqs. ij2U|) and ((2T|) be- 
cause of resonances around small rhombi. Let the addi- 
tional state with three parallel dimers on the parallelo- 
gram C2 be denoted by and the two dimer states on 
the perimeter of C2 be \ip2) and |03), ("02,3 l^i) = c, while 
("02 1 01 ) = c'. Then, the matrix elements E12 — E13 de- 
scribe short-loop resonances, while E23 corresponds to a 
tunneling along the composite loop. Generalizing trans- 
formation 1(19(1 to three states we finally obtain: 

Eu^-cEu- (d -^c 2 )E 23 , 

E12 ~ E12 , E23 ~ E23 — cEu (24) 

in the relevant case |c'| <C \c\ <C 1. Tunneling of dimers 
along a longer loop E23 is renormalized by short-loop 
hopping. Using the above expressions to calculate reso- 
nance of singlet bonds on C2 and C3 one obtains that in 
both cases 

t2 = h = -E23 = ■ (25) 

The potential energy given by the second term in En in 
((24(1 is again extremely small V^/^l ~ 0.02 and can be 
completely neglected^ Further extension of the above 
calculations to longer loops show that tunneling matrix 
elements of four-dimer moves are rather small ~ 0.07ii 
and should also be neglected. 
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C. Second-order mapping 

Analysis starts again with calculation of the mean-field 
contribution from 7i 2 to the ground state energy (diago- 
nal matrix elements) for an arbitrary dimer state. Every 
dimer has a finite energy contribution from 7^2^ \J 2 
while all nondimer bonds receive mean-field contributions 
from7T!2a: — j^J 2 - The total mean-field energy does not, 
therefore, depend on a chosen dimer covering of a trian- 
gular lattice and is equal to 




where N is number of sites on a kagome lattice. The 
mean-field energy l|26[l is subtracted in the following from 
T~L\ + H2, such that all diagonal matrix elements vanish. 

The off-diagonal matrix element of Ti.2 for the shortest 
loop C\ has nonzero contributions from 11.2a'- ~~igJ 2 > and 
from 7i2c- — jg J 2 - Combining them with the first-order 
result Eq. (|23|) we obtain 

ti = -E 12 = -p(l- 7 -j) . (27) 
Similar calculation for the loops C2 and C3 yields 

Taking into account Eq. (|20|l the tunneling matrix ele- 
ment between orthogonal dimer states along the loop C2 
(C3) becomes 

15 / 1 \ 
t 2 = t 3 = -E 23 = --^j(l--j) . (29) 

Loop C4 also acquires a finite tunneling rate between the 
two dimer states in the second-order. The corresponding 
matrix element is, however, small E-12 0.003 J 2 . 

IV. DISCUSSION 

Signs of the tunneling matrix elements calculated in 
the previous section have certain arbitrariness^* The 
negative sign of the resonance matrix element for the 
shortest loop C\ can be changed to positive by a gauge 
transformation. For this, real singlet wave- functions have 
to be multiplied by complex factors i n <-+ n i,e-ni,o j where 
for a given dimer configuration n r counts the number of 
dimers on links pointing upwards and right, n^ e (ni t0 ) 
counts the number of dimers on links pointing upwards 
and left from sites with even (odd) vertical coordinates. 
Dimers on strictly horizontal bonds do not contribute 
to the phase factor. By this operation resonance moves 
along every C± loop pick up an extra (—1) factor, chang- 
ing t\ — » — 1\. At the same time, amplitudes of dimer 
tunneling along C2 and C3 loops do not change sign by 
the above gauge transformation: t 2 ,s — > £2,3- An effective 



QD Hamiltonian for the trimerized kagome antiferromag- 
net is, therefore, dominated by the kinetic energy terms 
for resonance moves between orthogonal dimer configura- 
tions \(fi Cn ) and \tp' c ) for every loop C„ of three different 
types n = 1-3 on an effective triangular lattice: 

Wqd=^-^cJ(<J , (30) 

in 

h -a , ( 1 -54 " — —^-h'Y 

Amplitudes for three-dimer tunneling processes have no 
significant smallness compared to the strongest resonance 
move: t 2 ^/t 1 w —0.31 for J< 1. The kinetic coefficients 
t n are differently renormalized by the second-order pro- 
cesses such that importance of three-dimer moves is fur- 
ther increased towards the isotropic limit: £2,3 Ai ~ —0.5 
for J = 0.5 and £2,3/^1 ~ — 1 as J — > 1. 

The QDM (|30|l with only two-dimer resonances has 
been studied via mapping to a frustrated Ising model in 
transverse field: 5 '- The ground state of this dimer model 
is believed to be a crystalline \/l2x VT2 state, which con- 
sists of locally resonating dimer pairs and breaks trans- 
lational symmetry of the lattice. Such a state should 
have a fully gapped excitation spectrum. Properties of 
the QDM l|3U|) with several competing dimer resonances 
have not been studied so far. Dimer resonances along 
longer loops, though not very small, frustrate each other. 
The ground state of the realistic model 1(50) 1 should not 
be, therefore, very far from the idealized model with t\ 
terms only. In particular, we expect that the ground state 
breaks certain lattice symmetries. The excitation spec- 
trum is also expected to be gapped unless a fine tuning 
of t n drives the system towards an Ising-type transition 
point between two crystalline states. 

In conclusion, the presented derivation of the QDM for 
a realistic Heisenberg spin model on trimerized kagome 
lattice illustrates a generation of small energy scales in 
frustrated quantum magnets. The dimer resonance ma- 
trix elements in (|30|l are given by small fractions of a 
weaker exchange constant, e.g., t\ w 0.047 J. In a wide 
temperature interval t\ <C T <C J the quantum spin 
systems is described by an RVB liquid of singlet pairs. 
At very low temperatures T < t\ a valence bond crys- 
tal probably replaces an RVB state. The gap between 
the ground state and the first excited singlet levels is a 
fraction of t\. It is extremely difficult to resolve such a 
tiny energy scale in the exact numerical diagonalization 
of small clusters. Similarly, the low-temperature regime 
T < <i might be beyond experimental reach for possible 
realizations of the spin-1/2 trimerized kagome antiferro- 
magnet. The dimer crystallization at T — is driven by 
local resonances. Therefore, variational mean-field type 
approachesi 2 *^ are not capable to describe the precise 
nature of the corresponding ground states. 
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